Astronomy & Astrophysics manuscript no. h4350 


February 2, 2008 


(DOI: will be inserted by hand later) 





The Globular Cluster System of NGC4374 



M. Gomez 1 and T. Richtler 2 



m 
O 
o 

(N 

> 
O 



> 
00 

oo 



CO 

o 

S3 

Pl,: 

6 

CO 



X 



1 Depto. Astronomfa y Astroffsica, P. Universidad Catolica de Chile, Casilla 306, Santiago 22, Chile 
e-mail: mgomez@astro . puc . cl 

2 Grupo de Astronomfa, Departamento de Ffsica, Universidad de Concepcion, Casilla 160-C, Concepcion, Chile 
e-mail: tom@coma. cfm.udec . cl 

Received / Accepted 

Abstract. We study the globular cluster system (GCS) of the giant elliptical NGC 4374 (M84) in the Virgo cluster using B 
and R photometry. The colour distribution is bimodal with peaks at B-R=l.ll and B-R=1.36, fitting well to those found in 
other early-type galaxies. The radial profile of the cluster number density is flatter than the galaxy light. Using the luminosity 
function we derive a distance modulus of fi = 31.61 ± 0.2, which within the uncertainty agrees with the distance from surface 
brightness fluctuations. Blue and red clusters show similar radial concentrations and azimuthal distributions. The total number 
of clusters is N = 1775 ± 150, which together with our distance modulus leads to a specific frequency of Sn = 1.6 ± 0.3. This 
value is surprisingly low for a giant elliptical, but resembles the case of merger remnants like NGC 1316, where the low specific 
frequency is probably caused by the luminosity contribution of an intermediate-age population. A further common property 
is the high rate of type la supernovae which also may indicate the existence of a younger population. However, unlike in the 
case of NGC 1316, one cannot find any further evidence that NGC 4374 indeed hosts younger populations. The low specific 
frequency would also fit to a SO galaxy seen face-on. 
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1. Introduction 

The study of globular clusters in early type galaxies has reached 
a state, where surprises have become rare when only individual 
galaxies are studied. Progress in understanding the relation be- 
tween the morphology (i.e. the specific frequency of clusters, 
their spatial as well as their colour distribution) of a GCS and 
the host galaxy properties normally emerges from analyzing 
larger galaxy samples (e.g. Kundu & Whitmore 2001; Larsen 
et al. 1200 11 1 However, from time to time, one encounters galax- 
ies which exhibit some peculiarity in their GCS, which one 
would like to understand in a more general framework. 

For example, it is well known that "normal" elliptical 
galaxies have specific frequencies higher than about 3 (see 
Sect. 13.61 for the definition and Elmegreen 119991 for a review 
of specific frequencies). NGC 1316, the brightest galaxy in the 
Fornax cluster, has nevertheless a low specific frequency of 
only ~ 0.9 (Grillmair et al. 1999] Gomez et al. EDTJU . Since 
it is known as a merger remnant, it is tempting to seek an ex- 
planation not in the low number of globular clusters, but in 
the high luminosity due to the presence of an intermediate-age 
stellar population which formed in the merger a few Gyrs ago. 
Indeed, Goudfrooij et al. ( 1200 H . by means of spectroscopy of 
the brightest clusters, identified several intermediate-age glob- 



ular clusters among them. Also the apparently high SN la rate 
(NGC 1316 hosted SN 1980N and SN 1981D) could indicate 
a strong intermediate-age population, given that the progeni- 
tor population of SNe la is suspected to have an age of a few 
Gyrs (Yungelson et al. [T995l Yoshii et al. [T996l McMillan & 
Ciardullo fT996l ) 

The question arises whether the combination of a low spe- 
cific frequency and a high SN la rate is ubiquitous or whether 
more examples can be found pointing to the possibility to use 
the specific frequency as an indicator for the presence of a 
younger population. 

Early-type galaxies, which hosted more than one SN la, 
are rare. Besides NGC 1316, NGC 4374 in Virgo is one of 
the few ellipticals with two or even three SNe events: 1957B, 
19801 (also labeled "intergalactic" due to its location between 
NGC 4374 and NGC 4406) and 1991bg. 

The GCS of NGC 4374 is not well investigated. Ajhar et al. 
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( 1994) presented VRI photometry of globular clusters for 10 
galaxies in Virgo and Leo ellipticals, among them NGC 4374. 
Due to the small field, properties like the specific frequency or 
the spatial distribution could not be addressed. 

In the light of a possible intermediate-age population, it is 
worthwhile to investigate in detail the GCS of NGC 4374 with 
a larger field. 
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Fig. 1. The observed fields in the three nights. NGC 4374 is at 
the center. The bright elliptical at the left is NGC 4406 (M86). 

Table 1. The details of the observations. See Fig.[0for the ori- 
entation of the fields. 



night 


field 


B 


R 


19/3/99 


central 


30, 120, 5 x 600 s 


30, 120, 4 x 600 s 


20/3/99 


northern 


30, 120, 4 x 600 s 


30, 60, 4 x 600 s 


21/3/99 


southern 


120, 5 x 600 s 


60, 6 x 400 s 



the average of 4 neighboring pixels. After this, we combined 
typically 5 frames with exposure times ranging from 400 to 
600 seconds in both B and R. Cosmics were removed by a cr- 
clipping algorithm. The combined frames show a seeing of 1 .7" 
in B and 1.5" in R. 

A median filter was then applied to the combined frames to 
model and subtract the galaxy light. 

For the photometry, we used DAOPHOT under IRAF. 
Daofind was run on each of these frames with a detection 
threshold of 3 times the <x of the sky level. We then performed 
PSF photometry with allstar. The output lists were matched to 
leave only objects detected in both filters. 

For these objects we used the stellarity index of SExtractor 
(Bertin & Arnouts I1996> for distinguishing star-like objects 
from galaxies. The stellarity index ranges between 1.0 (star- 
like) and 0.0 (extended). This classification, however, becomes 
progressively more uncertain for fainter objects. 

As we do not resolve globular clusters at the Virgo distance, 
no light loss due to a different shape of clusters and stars is ob- 
servable. An aperture correction was applied in the usual way 
to calibrate our PSF-photometry with that of the standard stars, 
which were measured using a much larger aperture. 

2.2. The calibration of the photometry 

Only the first of 3 nights was photometric. About 40 Landolt 
( 1992) standard stars were observed at airmasses from 1.0 to 
1.8, giving a total of ~ 100 datapoints. Their colour range ex- 
ceeds that of the globular clusters. We then performed aperture 
photometry with radii from 4 to 30 pixels. A curve-of-growth 
was constructed for each standard, and only stars whose instru- 
mental magnitudes converged with increasing aperture radius 
were selected for the calibration. 

The resulting transformation equations are: 



2. Observations, reductions and photometry 

2.1. Observations 

The observations were carried out at the 3.5m telescope 
at Calar Alto, Spain, run by the Max-Planck Institute for 
Astronomy, Heidelberg. The observation period was 19 to 21 
March, 1999. The instrument was the focal reducer MOSCA 
(www.caha.es/CAHA/Instruments/MOSCA/index.html) 
equipped with a Loral 2K x 2K CCD. The pixel scale was 
0.513 "/pix and the usable unvignetted field of view ~ 13' x 13' 
(1.5K x 1.5K). The filters in use were Johnson B and R. In 
the first night several frames, centered on NGC 4374, were 
acquired, as well as the Landolt fields SA98, SA101 and 
SA107 at different airmasses. The observation log is given in 
TableG] 

Two additional fields were also observed during the second 
and third night, with an offset of ~ 8' to the north and south, 
respectively, as shown in Fig. [2 

Bias subtraction and flat-fielding were performed using 
standard IRAF procedures, resulting in a flat-field accuracy of 
about 1%. The images were aligned using the centers of many 
bright stars as reference points. Bad pixels were replaced by 



b - B + 0.874( ±0.026) + 0.223(+0.021) • X 

-0.176(±0.021)-(B-i?) (1) 
r = R + 0.027(+0.022) + 0.096(+0.016) • X 

+0.046(±0.020) (B-R) (2) 

with b, r, B, R the instrumental and standard magnitudes in the 
B and R filter. X is the airmass. The rms of the fit was 0.029 
and 0.027 for B and R, respectively. 

The calibration has been compared with photoelectric aper- 
ture photometry of NGC 4374 published by Poulain ( 1988 1, us- 
ing five aperture sizes. The mean difference was ~0.02 mag for 
each band, without any systematic trend. Local standard stars 
were then defined in the field to calibrate the non-photometric 
nights. 

We also compared the R magnitudes of the cluster can- 
didates with about 100 clusters in common with Ajhar et al. 
(fT994b . The mean difference is 0.006, with a <x of 0.085. No 
trend with magnitude is apparent (Fig.[5J- Regarding colours, 
no straightforward comparison is possible. However, we can 
plot V-I colours of Ajhar et al. versus our B-R colours and ask 
whether this relation matches the (B-R)-( V-I) relation for galac- 
tic globular clusters. This is done in Fig. |5] The agreement is 
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Fig. 2. Comparison of common objects in the photometry of 
Ajhar et al. ( 1994 1 and the present work in the R-band. AR 
means the magnitudes quoted by Ajhar et al. minus ours. The 
mean difference is 0.006 mag with a <x of 0.085 mag. 
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Fig. 3. Colour comparison of common objects in the VRI pho- 
tometry of Ajhar et al. ( 1994 1 and our BR photometry (filled 
symbols). We chose the V-I colour because of its widespread 
use in the photometry of globular cluster systems. Overplotted 
are galactic globular clusters, where the dereddened colours are 
taken from the compilation of Harris <1996t . There is no obvi- 
ous offset between the clusters of NGC 4374 and galactic glob- 
ular clusters, demonstrating the quality of our B-R colours. 



quite satisfactory. A few outliers may be caused by photometric 
errors. 



iii) uncertainty(B), uncertainty(B-R) < 0.3 

iv) stellarity index > 0.35 

These criteria arise from the assumption that clus- 
ters around NGC 4374 should resemble the Galactic ones. 
Regarding the color range, we have used the McMaster data 
(Harris 1996 1 and chosen galactic clusters with reddening 
smaller than Eb-v = 1-0, resulting in the above colour limits. 

The colours were de-reddened and the magnitudes 
extinction-corrected using the maps of Schlegel et al. (1998 1. 
We adopt E B - V = 0.04, A B = 0.173 and A R = 0.107. In the 
following, only reddening corrected values appear. 

2.4. Completeness correction 

For the later determination of the globular cluster luminosity 
function (GCLF) we need to know what fraction of globular 
clusters is found in each magnitude interval. A common way 
to evaluate the completeness is to perform artificial star ex- 
periments. Many artificial "star-like" objects are added to the 
frames and the detection, photometry and selection procedures 
are performed in exactly the same way as with the real cluster 
candidates. The number of recovered and selected objects per 
magnitude bin, divided by the initial number of artifical clusters 
gives an estimation of the completeness correction. 

To get sufficiently good statistics, one needs a large number 
of artificial objects, but care must be taken not to increase the 
crowding. Therefore, consecutive tests with a small number of 
artificial objects are preferred. 

Using the known PSFs, we added 200 artificial clusters in 
steps of 0.1 mag, starting from B=20.0 down to B=26.0, and 
applied the same photometric treatment and selection criteria 
as we did for the cluster candidates. 

We expect the completeness to vary with galactocentric 
distance because of the radially dependent galaxy light back- 
ground. However, a decrease in the completeness was only 
noticed for clusters located between radii of 50 and 150 pix- 
els (25'.'65 and 76'.'95). At larger radii, the completeness cor- 
rections stays constant with increasing distance. Fig. 0] (upper 
panel) gives the resulting completeness in dependence on B- 
magnitude for a colour of B-R = 1.2. 

We estimated the dependence of the completeness on 
colour using clusters of colours B-R=0.9, B-R=1.2 and B- 
R=1.5. In all, 140000 clusters were added for each colour. 
Fig.0](lower panel) shows that the completeness for red clus- 
ters is about 0.3 mag fainter than for blue ones. This difference 
must be taken into account for the calculation of the radial pro- 
file of red and blue cluster candidates (Sect. EOV 



2.3. Selection of cluster candidates 

886 objects have simultaneously been detected in the B and R 
bands. However, a fraction of them might be foreground stars 
and background galaxies. To statistically select globular clus- 
ters among them, we have used the following set of criteria: 

i) B > 20.5 

ii) 0.7 < B-R < 1.8 



2.5. Background correction 

Despite the applied selection criteria, there remain point 
sources, which may be foreground stars or unresolved back- 
ground galaxies. The surface density profile levels out at a ra- 
dial distance of about 500"(see Fig. II It . We therefore have 
used the outermost parts of the two fields observed during the 
second and third nights (see Fig.Q- That is, two rectangular ar- 
eas of 12'8 x 4'3 at the top of the upper field and at the bottom 
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Fig. 4. Top: The completeness for three annuli centred on 
NGC 4374, as function of the B magnitude, for clusters hav- 
ing B-R=1.2. The annuli are described by the following inner 
and outer radii (in pixels): 50<r<150, 150<r<300 and r>300. 
Only the inner annulus has a noticeable lower completeness. 
Bottom: The mean completeness (averaged over the three an- 
nuli), for clusters having colours of B-R=0.9, 1.2 and 1.5. 
These are represented by the three solid lines from the left to 
the right. The shift with colour is evident. The two dashed lines 
indicate the magnitude at which the mean completeness falls 
to 0.5 for blue (left) and red (right) clusters. The difference is 
~0.3 mag. 

of the lower field are our background fields. The photometry 
for these fields was done in the same way as for the central 
field and it was calibrated using the local standards defined in 
the overlapping regions. 

Completeness tests were also run on these two fields. The 
two background populations were averaged concerning colour 
distributions, luminosity function and surface density. 

3. Results 

We first present the colour-magnitude diagram and the colour 
distribution of those cluster candidates identified simultane- 
ously in B and R. After that we study the morphological prop- 
erties of blue and red candidates separately. In these subsec- 



tions we consider only sources brighter than B=23.8, where 
the mean completeness is ~ 60% (see Fig.^Ji. 

For deriving the GCLF (Sect. 13. 51 and the specific fre- 
quency (S n), we used only the R frame. The reason for doing 
so is that we reach more than half a magnitude deeper using 
the R frames alone, which might be crucial for the detection of 
the turn-over magnitude (hereafter TOM). We keep the crite- 
ria regarding magnitude range, photometric error and stellarity 
index. 

Completeness and background corrections for the R frames 
have been performed in the same way as for the combined sam- 
ple. 

3.1. The colour-magnitude diagram 

The B-R colour-magnitude diagram is shown in Fig. |5] In the 
upper panel, cluster candidates brighter than B=23.8 and hav- 
ing a radial distance less than 500"are plotted. The lower panel 
shows the same diagram for the background population in an 
area which is smaller by a factor 1.3 than that for the upper 
panel. It is apparent that objects bluer than B-R=0.9 are mainly 
foreground/background objects. This is in good agreement with 
the fact that B-R=0.9 is also a limit for the galactic clusters, as 
can be seen from Fig. [3] 

3.2. Colour distribution 

The colour histogram is given in Fig.|6] The upper panel shows 
the colour distribution of all sources having a radial distance 
less than 500"(bold line). The dotted line is the background 
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normalised to the same area. The lower panel shows the colour 
distribution for the background corrected sample (solid line). 
For comparison the B-R distribution of the galactic clusters is 
shown as well (dashed line), showing those clusters, for which 
B-R photometry is available and which have reddenings less 
than 1 .0 (8 1 clusters from the compilation of Harris 1996 1. The 
little peak at B-R=0.85 can be a residual from an inappropri- 
ate background subtraction. Without having further clues, we 
do not regard these objects to be globular clusters. The galac- 
tic distribution has the peak at B-R=l.l, qualitatively agree- 
ing with the blue peak of the NGC 4374 clusters. Although the 
comparison in the red regime is made difficult by the fact that 
many metal-rich galactic clusters with higher reddenings are 
omitted, it seems that the distribution of NGC 4374 clusters is 
more extended to the red and hence presumably to more metal- 
rich clusters. The appearance of the histogram does not depend 
much on the bin center which is demonstrated by the dotted 
curve. This is the sum of the contributions from all clusters, 
where each one is represented by a Gaussian at the correspond- 
ing colour. The only parameter is the <x of the Gaussian which 
we choose to be 0.05 to match the bin size. 

There are not many B-R photometries of GCSs available, 
so a comparison only with 3 galaxies, NGC 1380 (Kissler-Patig 
et aL lTWTl . NGC 1 1 99 and NGC 6868 (da Rocha et al. l2002b is 
principally possible. Unfortunately, in all three cases the colour 
distribution was not background corrected, which may hamper 
the comparison. However, the blue peak in these colour distri- 
bution was found at B-R= 1 . 1 in all cases, resembling the galac- 
tic system and, as we will see, that of NGC 4374. 

Ajhar et al. i 1994b and Gebhardt & Kissler-Patig J1999I re- 
port V-I photometry for GCs in NGC 4374. Ajhar et al. pointed 
out a strikingly narrow V-I colour distribution in comparison 
with other Virgo and Leo ellipticals, but again a more detailed 
comparison is made difficult by their small number statistics 
(moreover, no background subtraction has been attempted). 
Gebhardt & Kissler-Patig analyse the first 4 moments in the V-I 
colour distribution of a sample of 50 GCSs. NGC 4374 shows 
one of the highest skewness parameters in their sample, i.e. its 
colour distribution is strongly weighted to the red. 

The general appearance of the colour distribution looks bi- 
modal. A bimodal colour distribution has been found for many 
GCSs (e.g. Ashman & Zepf fT992l Larsen et al. l200Tl Kundu 
& Whitmore 2001 1, often interpreted as a signature of two 
episodes of cluster formation, sometimes in the context of a 
merger scenario. However, see Dirsch et al. J2003) for the ef- 
fects of a non-linear colour-metallicity relation. 

Following a common statistical approach, we performed a 
KMM-test (Ashman, Bird & Zepf ll994t on the colour distri- 
bution. This test is applied to the distribution uncorrected for 
the background population, in both homoscedastic and het- 
eroscedastic versions. It returns that the colour distribution is 
best represented by two Gaussians with positions at B-R=l.l 1 
(the blue peak) and B-R=1.36 (the red peak), and with re- 
values of 0.08 and 0.13, respectively (heteroscedastic mode). 
In the homoscedastic case, the peaks are at 1.14 and 1.42 with 
<x = 0.10. The P-value, which gives the probability of hav- 
ing a unimodal distribution, is practically zero in both ver- 
sions, i.e. a unimodal distribution is excluded. Fig.0shows the 
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Fig. 6. Top: the B-R colour distribution of cluster candidates 
around NGC 4374. The bin size is 0.05 mag. The thick solid 
line represents the colours before subtraction of the background 
population, whose distribution is shown with the dotted line. 
Bottom: the B-R colour histogram after subtraction of the 
background population, together with a so-called "generalised 
histogram" (solid line). This is constructed by placing a gaus- 
sian in the abscissa at each B-R colour, and then adding the 
contribution from all gaussians together. The latter is indepen- 
dent of the bin center, and in our case agrees with the "nor- 
mal" histogram. The dashed line shows the comparison with 
the galactic clusters. 



colour distribution together with the returned Gaussians in the 
heteroscedastic fit. In all GCSs investigated so far, the peaks 
exhibit approximately constant colours. Weak dependence on 
galaxy luminosity has been found for the blue and red peaks 
by Larsen et al. (2001 1, while Forbes & Forte (2001 ) see only 
a dependence in case of the read peak. 

As mentioned, also NGC 4374 fits to this notion compared 
with GCSs, for which B-R photometry is available. However, 
most work has been done in V-I (e.g. Larsen et al. 2001 Kundu 
& Whitmore 2001 ) where the peaks are located at V-I=0.95 and 
V-I=1.15. As can be seen from Fig. [3] these values fit well to 
the B-R peak colours quoted above. 

Is there a colour gradient? Fig.[8]plots the B-R colour ver- 
sus the projected distance from the galaxy centre (in logarith- 
mic scale). The scatter is large, but a slight trend that the clus- 
ters become bluer at larger radii, may be recognizable. 
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Fig. 7. This plot shows the colour distribution of the GCS of 
NGC 4374 (not background corrected) together with the two 
Gaussians, which best represent the distribution in the sense 
of a KMM-test. Only clusters redder than 0.925 have been in- 
cluded in the test. 
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Fig. 8. The colour gradient for cluster candidates brighter than 
B=23.8. The upper panel shows the B-R colour of individual 
clusters vs. the projected distance to the optical center. The B- 
R colour of the galaxy is shown for comparison (triangles). In 
the lower panel the mean color in concentric anuli (each 51'.' 3 
wide) is plotted. The error bars are the cr in the colour for each 
bin. The solid lines are the least-square fits. 



3.3. Blue and Red clusters 

We divided the clusters in two subpopulations at B-R=1.25 in 
order to search for possible morphological differences between 
the blue and the red clusters. 

The resulting counts are given in Tabled and the radial pro- 
files of red and blue clusters are shown in Fig. Rafter the sub- 
traction of the corresponding red and blue background levels. 
Red and blue clusters do not significantly differ in their radial 
distributions. They are well represented by: p(r) re( j x r~ L22±012 
andp(r) b i ue ocr- 100±015 . 
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Fig. 9. Radial profile of red (open circles) and blue clusters 
(filled circles), for objects brighter than B=23.8. They are not 
clearly distinguishable in their concentration. The solid lines 
are least-square fits. 

Table 2. The counts for the radial profile of blue and red sub- 
samples, using both B and R frames. The counts have been 
corrected for completeness and background contamination. 



r["] 


area [□'] 




yyrcd 


77/0 


6:890 


35.9 + 


5.8 


25.0 + 


5.1 


128'.' 3 


11.'483 


44.2 + 


6.3 


16.3 ± 


4.7 


179.' 6 


16: 076 


32.0 ± 


5.4 


20.4 ± 


4.9 


230/9 


20:669 


42.4 ± 


6.3 


19.7 ± 


4.9 


282/2 


25:262 


46.4 ± 


6.5 


15.4 ± 


4.6 


333'.'5 


29:856 


32.1 ± 


5.4 


17.1 ± 


4.8 



To investigate possible differences regarding the azimuthal 
distribution, we counted the cluster candidates in 12 equally- 
sized sectors, each spanning an angle of 30 degrees. The re- 
gion inside the radius r=51'.'3 (or, equivalently, ~ 5 kpc) was 
left out as the completeness is significantly lower. Clusters out- 
side r = 350" were not considered because the sectors become 
geometrically incomplete. The results (Fig. II Oi show that, on 
average, red and blue clusters are spherically distributed around 
NGC 4374. The galaxy's position angle of the major axis is in- 
dicated by the arrows at 135 and 315 degrees. 

There are some cases where blue and red clusters show sig- 
nificant differences in their azimuthal distributions, for example 
NGC 1316 (Gomez & Richtler l200TV NGC 31 15 (Kavelaars 
IT9981 . and NGC 1380 (Kissler-Patig et al. fT997l . In these 
galaxies, the red clusters resemble the shape of their host galax- 
ies, while the blue clusters are more spherical. The interpreta- 
tion as differences between halo and bulge/disk populations is 
suggestive. The elongations of the above galaxies, if caused by 
inclined disks, are strong enough to make a disk population 
of clusters easily distinguishable from a spherical halo popula- 
tion. This is not the case for NGC 4374 with its low ellipticity. 
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Fig. 10. The azimuthal distribution of the clusters. The posi- 
tion angle is counted north over east. The histogram at the top 
(solid line) represents the number of cluster candidates on each 
sector. The binning angle is 30 degrees. The dashed and dot- 
ted lines are the angular distribution of blue and red clusters, 
respectively. No systematic difference exists and both red and 
blue clusters show a circular distribution. The galaxy's posi- 
tion angle of the major axis is indicated by the arrows at 135 
and 315 degrees. 



Clearly, we cannot draw conclusions about halo and bulge/disk 
cluster subpopulations from their azimuthal distribution. 

3.4. Radial Profile 

We calculated the cluster radial profile by counting the clus- 
ter candidates in several annuli centered on NGC 4374, each 
25'.'65 wide, starting from a radius of r =8978 up to r =34674, 
and 10276 wide thereafter. The innermost ring (r =7770) was 
given a width of 5173 because of the low statistics. The results 
are listed in Table [3] The bin center and size are the first two 
columns. The area of the annuli is given in the third column. 
The fourth column shows the number of candidates up to the 
60% detection limit, before correcting for the completeness. 
The corrected number of clusters and its density are given in 
columns five and six. The remaining columns are used in the 
derivation of the specific frequency S n (see Sect. 13. 6> . 

The radial profile in R is shown in Fig.[^ The upper panel 
plots the GC number surface density and it is apparent that the 
cluster population extends at least out to r = 350". At larger 
radii the background dominates the profile. 

The background level is well determined and its value of 
8.36 ± 0.60 arcmin~ 2 has been subtracted from the tabulated 
densities, to derive the true radial profile, which is shown in the 
lower panel. The galaxy light in the R-band, arbitrarily shifted, 
is also plotted. The solid lines are: p(r) gcs oc r -i-09±0.i2 an( j 
p( r )gai K r _1,67±0 ' 02 , with p(r) gcs being the cluster surface den- 
sity (in units of arcmin~ 2 ) and p(r) ga i the surface brightness of 
the galaxy light. 

It is also apparent that the radial profile of the cluster den- 
sity around NGC 4374 is flatter than the galaxy light. 
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Fig. 11. The radial profile of the clusters' surface density for 
candidates brighter than R=23.5. Top: the number of candi- 
dates per square arcmin vs. the mean projected distance to the 
optical center. Bottom: the radial profile, after the correction 
for background objects. The galaxy light, arbitrarily shifted 
(crosses) and the least-squares fits (solid lines) are also plot- 
ted. 

3.5. The Luminosity Function 

The GCLF has been widely used as a distance indicator (see 
Hams l200Tl and Richtler 127)031 for recent reviews). GCLFs of 
a broad range of galaxy types are well characterised by log- 
normal or fs functions, the latter being of the form: 



t 5 (m) - 



3 VI: 



ncr, 



1 + 



(m — moY 
5oj~ 



(3) 



Ideally, the photometry should reach deeper than mo, the 
turnover magnitude (TOM). Unfortunately, our combined pho- 
tometry in B and R turns out to be insufficient to reach the ex- 
pected TOM for the Virgo cluster. The results of our artificial 
stars experiments show that we go half a magnitude fainter us- 
ing only R instead of both B and R. For this reason, we restrict 
to the R-band in the derivation of the GCLF. 

The counts are listed in Table [4] We remark that a gaus- 
sian fit does not give systematically different results (see e.g. 
Delia Valle et al. fT998l Larsen et al. l200Tl Kundu & Whitmore 
200l|, but a ts function leads to a better representation of our 
data. Furthermore, we assume a cr, of 1.1 which is a typical 
value for bright ellipticals (Harris 2001 Larsen et al. 2001 1. 
The least-square fit is shown in Fig.^1 where the data are the 
number of cluster candidates (corrected for completeness and 
background contamination) at each 0.5 magnitude in R. The 
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Table 3. This table shows the result of the radial profile for the surface density of cluster candidates. The first three columns list 
the center of each annulus (in arcsecs), its size and its area, in square arcmin. The observed number of cluster candidates brighter 
than R=23.5 on the corresponding annuli is given in column 4. After correcting for completeness, the actual number is shown 
in column 5. Column 6 lists this number divided by the area of the annulus, or the surface density of clusters. The next three 
columns are used for the derivation of the specific frequency S They are the number of candidates up to the observed TOM 
before and after the completeness correction (columns 7 and 8), and the total number of clusters in each annulus, assuming the 
LF to be symmetric around the TOM (column 9). The background level has been subtracted for these latter counts. For radii 
larger than 333V5, the annuli contain only background sources, so these four last rows are left blank. 
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Table 4. The counts for the GCLF in the R band. The first column is the bin center. The second column gives the observed number 
of candidates. The completeness factors (in column 3) lead to the completeness-corrected number of candidates in column 4. The 
observed background population and background completeness are listed in columns 6 and 7. The completeness correction for 
the background field is slightly smaller than for the central field due to better seeing. Finally, after subtraction of the background 
population (scaled to the area of the clusters), the total number of clusters is given in column 8. 
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error bars account for poisson error in the binning process and 
for the completeness correction. According to this fit, the TOM 
in the R band is 23.56+0.15. 

Leaving cr, a free parameter, we get for the TOM(R) = 
23.84+0.35 and cr, = 1 .22+011. We take the absolute TOM(R) 
from Delia Valle et al. (IT9981 who quote M R = -8.14 ± 
0.07 for the galactic system. Another absolute calibration 
comes from the Andromeda nebula. Barmby et al. ||2P01 ) give 



R=16.40+0.14 for the apparent TOM of M31. Together with a 
distance modulus of fi = 24.38 ± 0.07 (Freedman et al. 2001 1, 
one gets Mr = -7.98 ± 0.15. We average these two values and 
adopt M s = -8.06 ±0.15. 

This leads to distance moduli of 3 1 .6 1 ± 0.20 (keeping cr, = 
1.1) and 31.90 ± 0.38 (leaving cr, free) for NGC 4374. In the 
following, we consider only the TOM with cr, = 1.1. 
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Fig. 12. The GCLF in the R band. The solid line is the least- 
square fit of a fs function (equation 3) keeping cr = 1.1. The 
derived TOM is 23.56 + 0.15. If cr is left a free parameter, the 
resulting TOM is 23.84 + 0.35 and the function is given by the 
dashed curve. The vertical line at R~ 23.7 indicates the mag- 
nitude at which the average completeness is 50% using only 
R-band measurements. 

The SBF distance modulus is fi = 31.32 ± 0.11 (Tonry 
et al. 2001 1. Our distance, although marginally higher, is still in 
agreement within the uncertainties. A slightly higher distance 
would anyway be expected because we used zeropoints from 
the Galaxy and M31, which contain a larger fraction of blue 
clusters. Larsen et al. (1200 1> find that the V-TOMs of blue clus- 
ters are brighter than those of red clusters by an average of 0.37 
mag in their sample of 17 early-type galaxies. This difference 
should be somewhat smaller in R. Our photometry does not re- 
veal the TOM definitely. However, it only could be fainter. If 
this were the case the difference between the SBF distance and 
the TOM distance would grow larger. There are some examples 
where the disagreement between SBF and GCLF distance can 
be related to the existence of an intermediate-age population 
(Richtler 2003 1. Therefore, a deeper photometry would be of 
interest. 



3.6. Specific Frequency 

The specific frequency relates the number of clusters belonging 
to a given galaxy with its luminosity in units of My - -15. It 
was first introduced by Harris & van den Bergh ( 198TJ as: 



S N = N- 10' 



0.4-(Afv+15) 



(4) 



with N the total number of clusters and My the absolute mag- 
nitude of the galaxy in the V band. 

Once the turnover has been measured, My is normally de- 
rived via a calibration with the Galactic GCS. Neither My nor 
N are trivial quantities to measure. As it is nearly impossible to 
observe the total population of clusters, N must be determined 
through rough extrapolations. Typically, the number of candi- 
dates brighter than the TOM is obtained (or estimated, if the 
TOM is not reached by the observations), and then the counts 



are doubled, assuming the LF to be represented by a symmet- 
ric function (a ts in our case). Matters are further complicated 
if the observations do not reach the background level, in which 
case, a spatial extrapolation is required as well. 

The radial profile indicates that our population of candi- 
dates extends out to a projected radius of ~ 6 arcmin. This 
region is fully covered by our frames, so we do not need any 
extrapolation down to the background level. 

Again, we have used only the R frame to count the clusters 
in several annuli, in almost the same way as with the radial pro- 
file. The only difference is that we count now up to the TOM 
(R=23.56). Table[3]rists the counts for Sn- The bin center and 
size are given in the first two columns. The area of the annu- 
lus (in arcmin 2 ) is in the third column. The seventh column 
shows the observed number of candidates up to the TOM. The 
completeness-corrected number of clusters in each annulus is 
in column 8. Finally, assuming the LF to be symmetric around 
the TOM and after the correction for background objects (8.36 
arcmin -2 ), the total number of candidates on each annulus is 
given in column 9. 

We have left out the inner region of the galaxy (r <51'.'3), 
because the completeness decreases dramatically as a conse- 
quence of the high noise associated to the galaxy light. To esti- 
mate the contribution of this region to the total number of clus- 
ters, we have used the power-law describing the radial profile. 
This gives 65 candidates per arcmin 2 , or 166 candidates in this 
region down to the TOM. 

Adding the contribution of all annuli gives N = 1775 ±150 
candidates distributed over a total area of 112.53 arcmin 2 . 
The distance modulus derived in the previous section is fi — 
31.61 + 0.2 and the total apparent magnitude of NGC 4374 in 
the V band is 9.11 (de Vaucouleurs et al. I199U . Using the 
colour excess of Eb-v = 0.04 (Schlegel et al. 19981, the ab- 
solute magnitude is My = -22.62 ± 0.2; we therefore derive a 
specific frequency of Sn = 1.6 ± 0.3. This is a very low value 
for an early-type galaxy of this luminosity class. A possible 
explanation is suggested in the discussion. 



4. Discussion 

NGC 4374 is similar to NGC 1316 in two aspects which 
could both indicate the presence of an intermediate age- 
population: a high rate of supernova (SN) type la and a low 
specific frequency of globular clusters. SN type la are be- 
lieved to have their progenitors in populations older than 1 Gyr 



(e.g. Yungelson et al. [1995 Yoshii et al. 1996 McMillan & 
Ciardullo 1996 1 and a low Sn value of 1.6 might result from 
the higher galaxy luminosity produced by an intermediate-age 
population with respect to older populations. 

In the case of NGC 1316, this has been already proven by 
the identification of GCs with ages of about 3 Gyr (Goudfrooij 
et al. 20011. These clusters probably formed in a major merger 
event which was accompanied by a starburst. The very high star 
formation rate during these mergers seems to provide suitable 
conditions for the formation of massive stellar clusters (Larsen 
& Richtler l2000l 
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For NGC 4374, spectra of clusters are not yet available, and 
broad-band photometry remains inconclusive regarding the ex- 
istence of intermediate-age populations. 

Integrated spectroscopy as well gives no indication: 
Terlevich & Forbes ( 2002 1 have catalogued the ages and metal- 
licities of ~150 field and cluster galaxies. Using the SSP mod- 
els from Worthey & Ottaviani ( 1997 1 and the H/3 index mea- 
sured by Goudfrooij et al. (1999 1, they quote an age of 11.8 
Gyr for NGC 4374. The quoted ages for merger remnants like 
NGC 1316 and NGC 5018 (host to the la supernova 2002dj) 
are 3.4 Gyr and 1.5 Gyr, respectively. NGC 5018 also shares 
with NGC 4374 the characteristic of a low specific frequency 
(Hilker & Kissler [T996b . 

However, NGC 4374 is a polar-ring elliptical (Bettoni et al. 
120011 note also that among the 10 polar-ring ellipticals, listed 
in that work, 3 are la host galaxies). It has a dust lane, which 
extends out to about 1 kpc, and is a prominent radio source as 
well (see Hansen et al. 1985 for details on the morphology of 
the dust lane). These properties might be related to previous 
merger processes. 

An interesting note regarding la host galaxies and GCSs 
can be made from the paper of Gebhardt & Kissler-Patig 
( 1999 1. These authors analyse the V-I colour distribution of the 
GCs of a sample of early-type galaxies. Their "skewness" pa- 
rameter measures the asymmetry of the colour distribution with 
respect to the mean colour. The two GCSs which are skewed 
strongest towards red (i.e. metal-rich) clusters both belong to 
la host galaxies (NGC 4536, NGC 4374) as well as does the 
fourth in this sequence (NGC 31 15). 

Since SNe la in "normal" elliptical galaxies are so rare, one 
might ask whether NGC 4374 is in fact a SO galaxy seen face 
on. This would alleviate the problem of the low specific fre- 
quency, and given the uncertainties, a specific frequency of 1.9 
(the upper limit) is not unusual for a SO galaxy (Harris [T99 II 
Kundu & Whitmore l200TV 

Looking for more peculiarities, it is striking that the ellip- 
ticity of the isophotes of NGC 4374 decreases with increasing 
radial distance rather than the opposite, which is commonly 
observed for ellipticals (Caon et al. . 1994 1. It is plausible that a 
face-on disk, which in projection in the inner parts is dominated 
by an elliptical bulge, may show this feature. 

In that case, one may expect more conclusive informa- 
tion from dynamical analyses. Gerhard et al. (2001 1 included 
NGC 4374 in their sample of 21 round early-type galaxies for 
which they have derived the mass profiles out to 1 to 1.5 ef- 
fective radii. According to Magorrian & Ballantyne (2001}, a 
flattened structure when seen face-on, but dynamically treated 
as spherical, shows a spurious radial bias in its orbit structure. 
NGC 4374, however, does not show any striking anomaly re- 
garding its anisotropy parameters. On the other hand, the ra- 
dial extension of the analysis is limited to 100 arcsec of radial 
distance, so data which reach further out are of high interest. 
Investigating the kinematics of globular clusters is the best way 
to pursue this suspicion. 



5. Conclusions 

We analysed the GCS of the giant elliptical NGC 4374, host of 
two or even three SNe la. This may lead to the suspicion that an 
intermediate-age population may be present in NGC 4374, and 
perhaps manifest in the GCS properties. However, the colour 
distribution does not show anomalous features. It is best repre- 
sented by a bimodal distribution with peaks at B-R=l.ll and 
B-R=1.36, resembling the colour distribution found in other 
early-type galaxies. A unimodal distribution is ruled out by a 
KMM test. 

The spatial distributions of the red and blue subpopulations 
do not show significant differences. 

However, the number of clusters in NGC 4374 turns out to 
be surprisingly lower than what is expected for its high lumi- 
nosity. This could indicate the presence of an intermediate-age 
population, as has been the case for other early-type galaxies 
with low specific frequencies. It is intriguing that these exam- 
ples (NGC 1316 and NGC 5018) also have been hosts of type 
la SNe. 

The other possibility is that NGC 4374 is a SO galaxy seen 
face-on, in which case the low specific frequency of 1.6 would 
fit to other objects of this class. One then would expect a disk- 
like kinematic behaviour of the metal-rich globular clusters or 
of a subsample . Measuring radial velocities of many globular 
clusters could provide further clues. 
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